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A non-perturbative theory is presented which allows to calculate the solvation free energy of po- 
larizable ions near a water-vapor and water-oil interfaces. The theory predicts that larger halogen 
anions are adsorbed at the interface, while the alkali metal cations are repelled from it. The den- 
sity profiles calculated theoretically are similar to the ones obtained using the molecular dynamics 
simulations with polarizable force fields. 
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There are a number of long standing mysteries in the 
fields of physical chemistry and biophysics. The Hofmeis- 
ter effect which has now been known for over 120 
years is, perhaps, one of the oldest and most puzzling 
ones. Hofmeister observed that different ions have very 
different effect on stability of protein solutions. While 
some electrolytes are very efficient at salting-out pro- 
teins, others lead to protein precipitation only at much 
larger concentrations. A related mystery, which is also 
very old, has to do with the surface tensions. Some hun- 
dred years ago Heydweiller [§| noted that adding a strong 
electrolyte to water leads to increase in the surface ten- 
sion of the water-air interface. While the dependence 
on the type of cation is weak, there is a strong varia- 
tion of the excess surface tension with the type of anion 
— the lighter halides lead to larger excess surface ten- 
sion than the heavier ones. The sequence is precisely 
the reverse of the Hofmeister one. Both effects are com- 
pletely unaccounted for by the current theories of elec- 
trolytes, which go back to the pioneering work of Debye 
and Hiickel (DH) 

Application of the DH theory to the study of interfacial 
properties of electrolyte solutions was initiated by Wag- 
ner and continued by Onsager and Samaras (OS) [aj. 
These approaches were based on the observation that the 
image charge induced at the air-water interface repels 
ions from the surface, creating a depletion layer. The 
Gibbs adsorption isotherm then leads to the conclusion 
that the surface tension of aqueous electrolytes must be 
higher than that of pure water. Unfortunately, there is no 
way to account for ionic specificity within these theories. 
Since the hydrated size of all halide ions is nearly the 
same — and this is the only parameter that enters into 
DH theory — the OS approach predicts that the surface 
excess should be independent of the type of ion. Re- 
cently there have been proposed some other approaches, 
but none have proven completely satisfactory [7|. 

Some clues to the failure of the DH and the OS theories 
started to appeared in the 1990s when the photoelectron 
emission experiments Q and molecular dynamics simula- 
tions with polarizable force fields showed that contrary to 
the common wisdom, there were ions present at the air- 



water interface [9|, llO| . The simulations found that while 
hard alkali metal ions such as Potassium and Sodium and 
small halides such as Fluoride [ill are repelled from the 
interface, the softer more polarizable anions such as Bro- 
mide and Iodine are actually attracted to it [l2|- Pres- 
ence of highly reactive halogens at the air- water interface 
of aerosol particles might_ help to explain the excessive 
rate of ozone depletion [l3| . 

In this paper a new class of electrolyte models will be 
introduced. Unlike the previous approaches, the polariz- 
ability of ions will be explicitly taken into account. The 
new theory is intrinsically non-perturbative — all the mo- 
ments of the ionic charge distribution, and not just the 
dipole, are taken into account. The calculated solvation 
free energies are used to obtain the interaction potential 
of polarizable ions with an interface and to calculate the 
anion and cation density profiles. 

Since the pioneering work of Debye and Hiickel, ions 
have been modeled as hard spheres with a point charge 
located at the center While perfectly reasonable for 
describing bulk properties of electrolytes, this approach 
is bound to fail when applied to polarizable ions near a 
dielectric interface. The reason for this is easily under- 
stood by considering the simplest model of a perfectly 
polarizable ion idealized as a conducting spherical shell 
with a mobile surface charge. When such an ion moves 
across a dielectric air (oil)-water interface, the charge on 
its surface shifts progressively from the exposed air/oil 
portion to the part that still remains hydrated. For per- 
fectly polarizable ions, energy cost of charge localization 
is very low and is easily compensated by the decrease in 
the cavitational energy [14] as the ion moves across the 
interface. 

To make the discussion quantitative, consider a polar- 
izable ion — modeled as a conducting sphere of radius a 
and charge q, see Fig. 1 — at an air (oil)-water interface. 
Both half-spaces will be treated as dielectric continuums 
with permittivities e^, and £<, for water and air (oil), re- 
spectively. To gain insight into the problem we first con- 
sider an ion with one of its hemispheres submerged in 
water and the other exposed to air, h — a, Fig. 1. This 
problem can be solved exactly, yielding a purely radial 
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FIG. 1: Ion of radius a at an interface 



electric field and the electrostatic self energy of the ion 
given by 



(1) 



We can also calculate the ratio of the surface charge on 
the two hemispheres, obtaining Qo/qw — ^oj^w This 
means that for an air- water interface with ^o/^w — 1/80, 
the fraction of the ionic charge located in air is only 1% of 
the total charge! For a perfectly polarizable ion half ex- 
posed to air, 99% of its charge remains hydrated! This is 
very different from non-polarizable ions, which under the 
same conditions will have half of their charge exposed to 
the low dielectric environment, at a huge electrostatic self 
energy cost. This is the reason why non-polarizable ions 
can not penetrate into the interfacial region. Evidently 
this is not the case for polarizable ions which can adjust 
their charge distribution to minimize the electrostatic self 
energy cost. Unfortunately, once we leave the symmetric 
case of an ion located half way across the interface, no 
exact solution is possible and approximate methods must 
be used. 

The self energy of a perfectly polarizable ion with its 
lower extreme situated a distance h from the interface, 
see Fig. 1, can be written as. 
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where the capacitance C = af[eo/^w, h/a) and f(x,y) is 
a scaling function. In the limit <C e^,, we can expand 
/ in powers of eo/^w, 
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The value of the scaling function f{0,h/a) determines 
the capacitance in the limit of vanishing air permittiv- 
ity. In this limit the electric field lines originating on the 
charge inside water must be tangential to the interface, 
so that the normal component of the electric field van- 
ishes. Even this problem, however, is difficult to solve 
analytically. Exact solution is possible, however, when 
h/a ^ 1. In this limit all charge is confined to a small 
spherical cap located inside water. The curvature effects 
can be neglected, and the problem reduces to finding the 



solution of a mixed boundary value problem in cylindri- 
cal coordinates: V^0(z,p) = 0, 4''{0,p) = for p > p™, 
and <?ii(0, p) = y for p < p^, where prime refers to the 
derivative with respect to z, V is the electrostatic poten- 
tial of the spherical cap, and pm is its radius, see Fig. 1. 
Mixed boundary value problems are notoriously difficult 
to study. Fortunately, this particular one can be solved 
analytically using the Hankel transform techniques [3]. 
We find 
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The capacitance of the spherical cap can now be calcu- 
lated to be Cc = Pm/TT. We note that this is just half the 
value of the capacitance of a disk of radius pm in vacuum. 
This result can be understood by considering a charged 
disk in front of a dielectric medium of very low permit- 
tivity. The image charge induced on the interface will 
then be of the same sign as on the disk, and in the limit 
fo ^ £w it will also be the same in magnitude. Thus, one 
needs only half the charge of the disk in vacuum to have 
the same potential. 

In view of the natural symmetry of the problem it is 
convenient to express Cc in terms of the angle variable 9, 
so that Cc = aO/ir and /(O, h/a) = 9/tt. Writing the ca- 
pacitance in terms of 9 extend the range of validity of Cc 
outside the limit of completely flat disc to larger spherical 
caps. In particular, for a particle located half way across 
the interface, h = a, 9 = 7r/2, we find /(0, 1) = 1/2, 
which agrees precisely with the exact result of Eq. ([T]). 
Furthermore, comparing Eqs. ^ and ^ with Eq. ([T]), we 
see that /i(l) = 1/2. We are now in position to write an 
approximate expression for the self energy of a perfectly 
polarizable ion moving across a dielectric interface. 



Usih) 



q 



1 



2e^a + 



(5) 



2e„ 



where 9{h) — i?e[arccos(l — ^)]. The real part of 
arccos(a;) is used in order to continue its validity into the 
regions h > 2a, where 9{h) — tt. For h > 2a the electro- 
static self energy reduces to Us ~ (7^/2e^a, which is the 
usual Born self energy of a bulk ion. In writing Eq. ([5]), 
we have approximated the scaling function fi{x) by a 
constant, fi{x) — 1/2. This is permissible, since when 
the ratio Eo/e-w 1, the prefactor of /i is very small and 
the precise value of fi{x) is not important — it is com- 
pletely dominated by the first term of Eq. ([3]). We should 
note, however, that although Eq. ([5]) is very accurate for 
7reo/(2e^) < 9 < 37r/4, and for h/a » 1, it does not de- 
scribe perfectly the crossover from the interfacial to the 
bulk regime. The reason for this is that Eq. ([5]) does not 
fully account for the image contribution to the electro- 
static energy at distances h > a. It is possible to include 
this corrections into the theory at the expense of more 
complicated expressions. In practice, however, we note 
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that the image contribution is screened very strongly [ij , 
with the characteristic length equal to half the Debye 
length U'i.m{z) « exp(-2z/^£))/(4e^z). Therefore, 
for concentrations of electrolyte above physiological ones 
150mM, the image contribution decays to zero after only 
a few Angstroms. For now, we shall, therefore, ignore the 
image effect in the crossover region. 

It is important to stress the fundamental difference be- 
tween Eq. ([5]) and a similar expression for non-polar hard 
ions. If a hard ion is located half way across the interface, 
it will have half of its charge exposed to the low dielectric 
environment. The electrostatic self energy cost for this is 
approximately ~ /{tga). This is almost two orders of 
magnitude larger than the energy cost for a polarizable 
ion to be at the same position! This is the reason why 
non-polarizable ions will not be found at the interface. 

The force that drives ions towards the interface arises 
from the cavitational energy. Presence of ions disturbs 
the hydrogen bond network of water and costs energy. 
We can estimate this energy cost by considering a cavity 
which must be formed in water to accommodate an ion. 
For small cavities of radius a < 4 A, which do not signifi- 
cantly perturb the hydrogen bonds, the energy cost scales 
with the volume of the void, while for larger cavities the 
energy cost scales with the cavity surface area [3l ■ This 
is, the so called, hydrophobic crossover from small to 
large length scales [17|. Small alkali metal and halogen 
ions are in the volumetric scaling regime. If one part of 
an ion leaves the aqueous environment, the cavitational 
energy will decrease proportionately to the volume ex- 
posed. This results in a short range cavitation potential 
which forces ions to move into the air (oil) phase. 
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From the results of bulk simulations we calculate 
that V w 0.3/cbT/A^. To account for the fact that a cav- 
ity containing an ion should be somewhat larger than the 
ionic crystallographic size, we will use v « 0.5fcsT/A^. 
This corresponds to the cavity radius about 20% larger 
than the crystallographic radius. 

For small non-polarizable ions, the cavitational energy 
is not sufficient to force ions into the low dielectric envi- 
ronment — the electrostatic energy cost is way too large. 
On the other hand, for soft polarizable ions, the elec- 
trostatic self energy cost is very small, since the ionic 
charge distribution can easily deform to remain mostly 
within the hydrated portion of the ion. The significant 
gain in the cavitational energy, and the low cost in elec- 
trostatic self energy, makes it energetically favorable for 
polarizable ions to move into the interfacial region. The 
amphiphilic nature of large ions, such as hexafluorophos- 
phate PFq , has been known for a long time. The cavita- 
tional energy for these ions is so large, that they actually 
adsorb to the interface, lowering its surface tension (18j . 



What has been discovered recently is that smaller po- 
larizable ions apparently can also have some amphiphilic 
activity — although not sufficiently large to lower the 
interfacial tension • 

So far we have considered only perfectly polarizable 
ions. Real ions, however, have finite polarizability. It is 
not obvious how the effects of finite polarizability can be 
included in the the formalism developed above. In fact 
it is not even clear, if the concept of bulk polarizability, 
as a liner response to the external field, is relevant in 
the interfacial geometry, where a rapid variation of the 
dielectric constant makes all the moments of the charge 
distribution — not just the dipole — relevant. For per- 
fectly polarizable ions we have avoided this difficulty by 
solving the full electrostatics problem for a conducting 
sphere. For ions of finite polarizability, to have a com- 
pletely quantitative picture it might be necessary to go 
to full ab initio calculations. In the absence of such, we 
can still gain some insight into this difficult problem by 
considering a simple model. In the spirit of Landau, we 
will construct the polarization energy Up by exploiting 
the symmetries of the problem. 

Consider an ion of radius a and bulk polarizability 7. 
We will define the relative polarizability of this ion as 



a = 7/70, where 70 



is the polarizability of an ideal 



ion of the same radius, modeled as a conducting sphere. 
For non- ideal ions with < a < 1, the surface charge can 
not fully adjust to the external electric field. Therefore, 
for such ions, there is an additional non-electrostatic — 
quantum mechanical — energy cost for dislocating ionic 
charge from its position of equilibrimn. Suppose that for 
a given ion the fraction of its total charge inside water is 
a;, then the charge exposed to air (oil) will be (1 — x)q. 
For ions which are highly polarizable x « 1, as long as 6 
is not too small . There is, however, a polarization energy 
cost for shifting a fraction of the ionic charge (assumed 
to be originally uniformly distributed along the surface of 
the ion) from its equilibrium position in the air portion 
of the ion to the water part. Within our simple dielectric 
model Up must be invariant under the transformation 
q —q. It must also be invariant under the transfor- 
mation 6 ^ TT — 9, when eo ^ and x 1 — x. To 
respect these symmetries, the polarization energy must 
be an even function of the difference between the initial 
(before exposure) and the final (after exposure) amount 
of charge on the part of the ion exposed to air/oil. Taking 
all these considerations into account and recalling Eq. ([5|) 
for a perfectly conducting sphere, we arrive at the polar- 
ization energy for a non-ideal ion. 



PUp{h;x) = ^ 
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where 13 = l/ksT and As = q^/ewksT is the Bjerrum 
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length in water. The terms in the first square brackets 
of Eq. ([6]) are the electrostatic self energy costs of the 
parts of the ion exposed to the water and the air, re- 
spectively. The second brackets contain the energy cost 
arising from the induced inhomogeneity in the ionic sur- 
face charge distribution. The coupling constant g must 
be a function of the relative polarizability, g{a). In the 
limit a — > 1, ion becomes perfectly polarizable, so that 
g ^ 0; while in the limit a 0, ion becomes completely 
hard and g oo. To account for these, we will write 
g{a) = x(l — ct)/ct, where x is a pure number. The pre- 
cise value of X can only be obtained from the ab initio 
calculations. For now we will take it to be of order unity, 
X w 1. Minimizing Eq. ([6]) with respect to x gives the 
fraction of the total charge located on the hydrated part 
of the ion, Xmin{h)- Substituting this back into Eq. 
yields the polarization potential that an ion feels as it 
moves across the interface Up{h) = Up{h; Xminih)). For 
an ideal ion of a = 1 located half way across the interface 
h = a, the energy Up(a) reduces precisely to the expres- 
sion given by Eq. ([T|). Therefore, for this case, the for- 
malism developed above gives the exact result. The total 
solvation potential felt by an ion of relative polarizability 
a is Usoi{h;a) — Ucav{h) + Up{h). Using this potential, 
we can calculate the ionic density profiles inside a small 
water droplet of radius R by explicitly solving a modified 
Poisson-Boltzmann equation, 

V (pir) = — p 

-/3g0(r)-/3(7,„,(/i;Q,) \ 
-R^ ' (7) 

where r is the distance measured from the center of the 
water droplet, h = R — r, and ac and aa are the relative 
polarizabilities of cations and anions, respectively. The 
density profiles for electrolyte solutions of KI, KBr, and 
KCl are presented in Fig. 2. The polarizabilities of ions 
7if = 0.79A3, 7j = 7.4A3, jbt = 5.07A^, jci = 3.77Af , 
where take from Ref. [3] and the ionic sizes ax — 1.49A, 
a/ = 2.05A, GBr = 1-8A, aci = 1.64A, from Ref. [26\. In 
agreement with the polarizable force fields simulations, 
the theory predicts that Iodine is strongly adsorbed at 
the air water-interface. We also find that there is a sig- 
nificant concentration of Bromide, while Chloride, Potas- 
sium and Fluoride (not shown) are depleted from the in- 
terfacial region. The current theory, however, predicts 
that there is less halide ion adsorption than is found in 
the simulations. The difference might due to the over- 
estimate of the neat water-vapor surface potential pre- 
dicted by the polarizable force field simulations to be 
— 500mVR:! —20kBT/q. Such a huge junction potential 
will irreversibly drive polarizable halides toward the va- 
por phase, resulting in a large density build up along the 
interface. Recent ab initio simulations |21l] . however, find 



a much smaller contact potential, — ISmYw —0.7kBT/q, 
for a water-vapor interface. This might lead to a smaller 
adsorption, in line with the predictions of the present 
theory. 

There is a subtle interplay between ionic size and po- 
larizability. Although the absolute value of polarizability 
for halide ions varies dramatically, the relative polariz- 
ability a is practically the same for all ions ^ 0.85. On 
the other hand, the relative polarizability of Potassium 
is almost four times lower. Future work should try to 
optimize the parameters of the theory v and x by fitting 
Usoi to the full ab initio simulations. Unfortunately, no 
such calculations are available at this moment. Finally, 
since the surface tension of electrolyte solution is directly 
related to ionic adsorption in the interfacial region, the 
theory developed also accounts for the Hofmeister series 
for halogens. 
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anion directly at the interface results in a build up of K in its vicinity. 
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